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ABSTRACT 

SN 200 lem, identified as a Type Ic supernova, has recently been detected in the radio and X-rays, 
> 2 yr after the explosion. The high luminosities at such late times might arise from a relativistic 
jet viewed substantially off-axis that becomes visible only when it turns mildly relativistic and its 
emission is no longer strongly beamed away from us. Alternatively, the emission might originate from 
the interaction of the SN shell with the circumstcllar medium. We find that the latter scenario is 
hard to reconcile with the observed rapid rise in the radio flux and optically thin spectrum, F v oc 
j,-o.36±o.i6£i.9±o.4^ w i)£ e th e s e features arise naturally from a misaligned relativistic jet. The high 
X-ray luminosity provides an independent and more robust constraint - it requires ~ 10 51 erg in 
mildly relativistic ejecta. The source should therefore currently have a large angular size (~ 2 mas) 
which could be resolved in the radio with VLBA. It is also expected to be bipolar and is thus likely 
to exhibit a large degree of linear polarization (~ 10% — 20%). The presence of a relativistic outflow 
in SN 2001em would have interesting implications. It would suggest that several percent of SNe Ib/c 
produce mildly relativistic jets, with an initial Lorentz factor Tq > 2, while the fraction that produce 
GRB jets (with Tq > 100) is ~ 100 times smaller. This could considerably increase the expected 
number of transients similar to orphan GRB afterglows in the radio, and to a lesser extent in the 
optical and X-rays, if there is a continuous distribution in Tq. Furthermore, this may give further 
credence to the idea that core collapse SNe, and in particular SNe Type Ib/c, are triggered by bipolar 
jets. 

Subject headings: stars: supernovae - supernovae: individual (SN 2001em) — gamma-rays: bursts — 
ISM: jets and outflows 



1. INTRODUCTION 

Supernova (SN) 2001em was discovered on Septem- 
ber 15, 2001 in the nearby galaxy UGC 11794 
ijPapenkova et alJl200l(l . at a redshift of z = 0.019493. 
This corresponds to a distance of D w 80 Mpc (for 
Oa = 0.7, Qm = 0.3 and h = 0. 71). It was classified as a 
Type Ib/c SN (most likely Ic, iFilippenko fc Chornockl 
120011) . SNe Type Ib/c - some of which are thought 
to arise from the core collapse of a Wolf-Rayet (WR) 
star - have drawn more attention in recent years due to 
their association with gamma-ray bursts (GRBs). The 
best and most secure as sociation so far is between GRB 
030329 and SN 2003dh (jStanek et al.ll200l IrTiorth et all 
2003). A compelling case also exists for SN 1998b w (at 
Z = 0.0085) and GRB 980425 (|Oalama et alJli^) . 

This raised interest in the search for signatures of GRB 
jets in nearby Type Ib/c SNe (e.g.. iPaczvhskil [2001). 
Typically, the narrow GRB jets point away from us and 
will not be detectable in 7-rays, but the SN might still 
be observed. As the off-axis GRB jets become mildly 
relativistic, months to years after the explosion, their ra- 
diation is no longer strongly beamed away from us, and 
they could become detectable in the r adio. 

Thus motivated, iStockdale et alJ l|2004|) observed a 
large sample of SNe Ib/c at late times, and detected 
SN 2001em on October 17.18, 2003 at 8.4 GHz as a 
1.151 ± 0.051 mJy radio source. In addition to its high 
radio luminosit y, Lr ~ 10 28 erg s" 1 Hz" 1 (second only 
to SN 1998bw; IKulkarni et aTlri998|> . SN 2001em was 
also unusual in its subsequent evolution. Its 8.4 GHz 
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flux rapidly increased to 1.480 ± 0.052 mJy on January 
30.90, 2004. This corresponds to a temporal index of 
a = 1.9 ± 0.4, where F v oc v^t a . Interestingly, the 
source appeared nonthermal, exhibiting a spectral slope 
of (3 = -0.36 ± 0.16 between 4.9 and 14.9 GHz, at the 
second epoch. More recently, on April 4.81 2004, Chan- 
dra detected SN 2001em in the X-ray (0.5 - 8 keV) with 
a luminosity of Lx ~ 10 41 erg s _1 and j3 ~ —0.1 ± 0.35 
llPoolev et al.H 200l. 

In this Letter we investigate different explanations for 
the unusual emission from SN 2001em. The two most 
natural mechanisms are (i) the interaction between the 
SN shell and the circumstellar medium (CSM), and (ii) 
off-axis relativistic jets. We examine these two possibili- 
ties in detail in SjHand respectively. Our conclusions 
are discussed in 21 

2. INTERACTION BETWEEN THE SUPERNOVA SHELL 
AND THE CIRCUMSTELLAR MEDIUM 

The characteristic SN radio lightcurves are thought to 
arise from the competing effects of a slowly declining 
nonthermal radio emission and a more rapidly declining 
absorption. Under the assumption that the fractions of 
internal energy in magnetic fields (es) and in relativistic 
electrons (e e ) remain constant with time, the observed 
radio flux can, to first approximation, be written as F v oc 
vPt a e~ T ', where (3 = — ^ and p is the power-law index of 
the electron energy spectrum ()Chevalierll994|L The early 
optically thick phase, r > 1, can be dominated by either 
free-free absorption or synchrotron self-absorption. The 
high expansion velocities and low CSM densities found in 
type Ib/c SNe suggest that synchrotron self-absorption 
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is the dominant mechanism in these objects l)Chevalierl 
1998). Synchrotron self-absorption leads to a power law 
in both time and frequency, F v oc i/ 5 / 2 , instead of an 
exponential form for free-free absorption. 

SN 200 lem showed both a fast rise in its radio flux, and 
an optically thin spectral slope, F v oc Z y-°- 36 ± - 16 i 1 - 9 ±0- 4 . 
While the former may be similar to that expected from 
synchrotron self-absorption, the latter is clearly not. The 
fact that the rapid rise occurs with little absorption, im- 
plies that it is not because of a reduction in optical depth. 
The usual models described above therefore fail to re- 
produce the observed increase in flux. This behavior has 
not been observed previously in radio SNc, although SN 
1987A (Type II) has shown a strong rise in its radio flux 
(Ball et al.lll995j) toge ther with an optically thin spec- 
tral slope, f3 w -0.95 l|Manchester et al.l f2002L that has 
been attributed to interaction with the dense w ind from 
a previous evolutionary phase ijChev alicr 1993). 

In order to address the question of whether or not the 
radio emission seen in SN 2001em is consistent with syn- 
chrotron radiation from the interaction of the SN shell 
with th e CSM, we generalize the analysis of iWaxmanl 
l)2004a|) . which applies to expansion in a p cx t oc r~ 2 
medium, to p cxt = Ar~ k . Let us consider a sub- 
relativistic shell ejected by the SN explosion, with mass 
M, total energy E, and initial velocity vq. Denoting, 
£dcc, the time at which the SN shell decelerates sig- 
nificantly, we have t dcc = [2(3 - ^E/AttAv^} 1 ^ 3 -^ 
and v pa u x min[l, (t/t dcc ) ( - k ~ 3 ^ ( - 5 ~ k ^}. The sharp rise, 
a = 1.9±0.4, and the spectral slope, (3 = — 0.36±16, that 
were observed in SN 200 lem, cannot be achieved after 
idee (jFrail et al.ll2000f l. On the other hand, at t -C ^dec, 
the observed spectral slope suggests that we are in the 
power law segment of the spectrum where (3 = —z^-, 
which implies a = 3 — M5±£i. j n order to obtain a ps 1.9 
one needs k < 0.55 — 0.63 for 2 < p < 3. Such a smooth 
power law density profile is u nlikely in the immediate sur - 
roundings of a massive star ijGarcia-Segura et alj [l996). 

Explaining the X-ray luminosity, Lx ~ 10 41 erg s _1 
at t w 950 days, is not trivial. We have L x ~ 
fx£r a .d£e(E/t)mm{l 7 (t/t dcc ) 3 - k ] where f x is the frac- 
tion of the radiated energy in the 0.5 — 8 keV Chan- 
dra range, e ra d ~ min[l, (7 m /7 c ) p ~ 2 ] is the frac- 
tion of the energy in electrons that is radiated away, 
and min[l, (t/tdec) 3_fe ] is the fraction of the total en- 
ergy E that is in the shocked CSM. This implies 
(3/x)(10e rad )(3e e ) J B 5 i min[l, (^cc/10 3 d) fc - 3 ] ~ 1, where 
E 51 = E/(10 51 erg), which suggests that 1 E 51 > 1 
and idec S 10 3 days. The latter condition implies 
V c £ 0.5(^ 5 iM*) 1/3 for k = 2, where A* = A/(5 x 
10 11 gr cm -1 ). As a consequence, the velocity of the 
ejecta must be at least mildly relativistic with E51 ~ l. 2 

The extrapolated radio flux in 8.4 GHz at the time of 
the X-ray observation is ~ 1.7 mJy, which corresponds 
to a radio luminosity of Lji ~ 10 38 erg s _1 . This would 
lead to (3 ps —0.6 for a single power law in that en- 
ergy range, which is consistent with p ps 2.25, as long 

1 The bare minimum for the energy content is E ~ 10 49 erg 
for fx e rad e e — 1- Such an extreme efficiency is, however, highly 
unlikely. For more reasonable values of fx , £e ~ 1/3 and e ra[ i ~ 
0.1, we need E ~ 10 51 erg. 

2 In this case, only a small part of the mass in SN shell, M ~ 
E/c 2 ~ 5 X 10~ 4 E5i Mq, would have an initial velocity vg ~ c. 



as v c > 10 18 Hz. The ratio Lr/Lx requires p < 2.25, 
where p < 2.25 gives v c < 10 18 Hz and v c {p — 2) ~ 
10 16 Hz. Such high values of v c favor a low CSM density, 
A* < 0.03(3e B )- 1 (l+r)- 4 / 3 , where Y is the Comptony- 
paramctcr which satisfies Y(l + Y) ~ (v/c)e ra de e /es. In- 
terestingly, a similarly low value of A* is required in order 
to explain t he lack of detection of an off-axis GRB jet in 
SN 1998bw l!Waxmadl2004albT:ISoderberg et alJl2(¥)l . 3 

3. EMISSION FROM AN OFF-AXIS RELATIVISTIC JET 

We first consider the off-axis emission from a uni- 
form double sided jet with an initial half-open ing an- 
gle 9 and sharp edges (e.g.. IGranot et af]l2002|) . Later, 
we briefly address 'structured' jets, where the energy 
per solid angle, e, smoothly decreases with the angle 9 
from the jet symmetry a xis, e oc 9~ 2 (jRossi et al.ll2002l 
iZhang fc MeszaroJ20f)2ft. 

Following G ranot fc Loebl l)2003ft and generalizing their 
results to a stellar wind external density profile, p cx t = 
Ar~ 2 , wc obtain expressions for the radius Rj where the 
Lorentz factor 7 of the jet drops to Bq , and the radius 
-Rnr where the jet becomes sub-relativistic, 

Rj = Rnn/f = E/2ttAc 2 = 3.5 x 10 17 E^A' 1 cm , (1) 

where / « 1 - hi9 a , and E — lO 51 ^! erg is the en- 
ergy of the jets. The typical angular s ize of the jet at 
the n on-relativistic transition time (see iGranot fc Loebl 
EH, t NR ~ Rnk/c, is 

where Da is the angular distance to the source. At the 
distance of SN 2001em, NR = 0.88(//3)£ 5 i A7 1 mas. 

The temporal index a ~ 2 is consistent with the ris- 
ing part of the light curve for a GRB jet viewed off-axis 
from an angle of 9 n hs a few 9r> w.r.t. the jet axis 
(e.g., Fig. 2 of IGranot et~aT1 12002). Therefore, wc ex- 
pect the peak flux to occur at i pca k = 3C yr, where 
C > 1. The peak flux at v — 8.4 GHz should be about 
i^.peak ~ 2C 2 mJy. Given the late peak time, it is 
likely that 9 Q \, & > 1 and therefore i pca k ~ ^nr ~ 3 yr, 
for which the source angular size is ~ #nr- Accord- 
ing to Eqs. & |3 #nr ~ 2 — 3 mas. One can also 
estimate #nr by requiring an apparent velocity of c, 
#nr ~ ct^/DA ~ 2.4(^nr/3 yr) mas. Such an angu- 
lar size could be resolved by VLBA. 

In order to explain the spectral slope of (3 ~ —0.4, we 
require that v m < v < v c , for which f3 = — The 
measured value of (3 can be somewhat lar ger than this 
asymp totic value if v m ~ 1 GHz . Following iNakar et alJ 
Ipf)^ and IGranot fcS^rll ll200l . we find 

P -QQt: Si?) „- Pr p-l >+l)/4 ,3(p+l)/4 

■fu,peak — g(2 2) £ e,-l e S,-2 A * 

t pcak = a (^j tj = 34(1 + z)aj±el hs days , (4) 
(for 6» obs > 26» and at the redshift of SN 2001em) where 
g (p) = (p - 0.18)e" 166 P (|Ef) P » and a elates tj to 

3 IWaxmanl I2004H) also derived v ~ 0.8c for SN 1998bw, al- 
though with a relatively low energy of E ~ 10 49 erg. 
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ipcak- For ^obs ~ 1 we expect a ~ 4, while for # b s <C 1 
we expect a ~ 1. For SN 2001cm, F^peakClO GHz) w 
2C 2 mJy, which equals the flux in Eq. 01 for a w 4 
and o bs ~ C _1 / p (7r/2). This suggests a viewing angle 
fobs > C _1 / p rad. Since for SN 2001em we know that 
ipcak = 3C yr, Eq. H yields aE^A^O 2 ^ ~ 30C, which 
implies a small CSM density, A„ ~ 0.1, similarly to [2J 
This relation can also be used to simplify Eq. [3] and 
eliminate the dependence on 9 b s and a, 

F„ peak ~ o.2c-^r-i4 p+ - 2 )/4 ^i 3 " p)/4 4r l)/2 ^^ p)/2 mJy • 

(5) 

Thus we obtain that e^^A^^E^/ 2 ~ 
10C P . Assuming a typical energy in GRB jets of -E51 ~ 1, 
and A* ~ 0.1, this gives <?~\^H^ ~ 6C P . As dis- 
cussed in the ratio Lx/Lr implies p < 2.2. For 
C«l the above condition can be readily satisfied for a 
wide range of reasonable parameter values (e.g., e e ~ 0.3, 
€b ~ 0.1). However, since e e , es < 1/3 — 1/2, we must 
have C < 2 — 3, which implies i pea k ^ 5 — 8 yr. 

Finally, we briefly address a 'sructured' GRB jet 
viewed from a large angle # bs- If the jet has an outer 
edge at # max < 9 ohs , then the light curve would not be 
very diff erent from that fo r a uniform jet viewed at o b s > 
80 fe.g.. lWei fc Jinll2003|) . In this case, the above anal- 
ysis is still approximately valid. If, on the other hand, 
#max = tt/2 or o b s < #max, then the early light curve 
is dominated by emission from material along the line of 
sight. In this sharp rise like the one observed in 

SN 2001em (a « 2), together with the observed spectral 
slope, (3 ~ —0.4, cannot be achieved after the time i<jec 
when the material along the line of sight decelerates sig- 
nificantly iGranot fc Sarill2002HKumar fc Granotl 120031: 
IGranot fc Kumarl 12003ft . Therefore, the only way this 
scenario might still work is if we are before td cc - In this 
case <p Cak > 3 yr is given by t dcc ~ ^nrF^" 2(4 ~ fe)/(3 ~ fc) . 
This suggests a mildly relativistic initial Lorentz fac- 
tor along the line of sight, To < a few, which might 
also explain why no GRB or X-ray flash w as observed 
(e.g.. lRamirez-Ruiz fc Llovd-Ronning|l2002|) . despite the 
very low redshift of SN 2001cm. Similarly to the non- 
relativistic case discussed in ^ a ~ 2 requires k < 0.6, 
which is unlikely. 

4. DISCUSSION 

Different possible explanations for the radio emission 
from SN 2001cm > 2 yr after the SN have been con- 
sidered. We find that the large temporal index, a = 
1.9 ± 0.4, together with the optically thin spectral slope, 
(3 = —0.36 ±0.16, cannot be naturally explained as emis- 
sion from the interaction between the SN shell and the 
CSM. This would require either an a lmost uniform exter- 
nal d ensity, or a density bump fe.g.. lRamirez-Ruiz et abl 
120011) . On the other hand, we find that a GRB jet, or 
even a jet with a mildly relativistic initial Lorentz factor, 
To > 2, that points away from us can naturally reproduce 
the observed temporal and spectral properties. 

Since the actual observed rise in the radio luminosity 
was only ~ 30%, it might still not be indicative of a long 
episode of increasing flux and could be only due to a 
local density bump. However, the measured X-ray lumi- 
nosity provides a stronger and more robust constraint. It 
requires ~ 10 51 erg in ejecta with a mildly relativistic ex- 



pansion velocity. Such a system would be physically very 
similar to an initially relativistic jet which became mildly 
relativistic at t ~ £nr, and began to approach spherical 
symmetry. It is also reasonable to expect that the mildly 
relativistic SN ejecta would be somewhat elongated along 
the rotational axis, similar to a relativistic jet near £nr- 
A high degree of linear polarization might therefore be 
expected. The polarization from a relativistic jet viewed 
off-axis is expected to reach its maximum value near the 
time of the peak in the light curve, t pC ak- For a rela- 
tivistic jet the peak polarization can reach ~ 30% — 40%, 
while for a mildly relativistic jet it is probably more mod- 
est, ~ 10% — 20%, but still significantly higher than for 
a typical SN. 

The best way to test our conclusion of a mildly rela- 
tivistic expansion velocity is via the angular size of the 
image, which should be > 2 mas, and could be resolved 
with VLBA. For a double sided relativistic jet, we might 
observe both jets, if the viewing angle is large enough, 
$obs ^ lj so that the difference in brightness between th e 
two jets would not be very large l|Granot fc Loebl l2003). 
In this case, their brightness ratio and its temporal evo- 
lution can help determine our viewing angle, 9 a b s - 

If indeed the radio and X-ray emission observed in SN 
2001em are from an off-axis relativistic jet, then this 
has several interesting implications. This could provide 
an estimate for the fraction /rj of SNe Ib/c that pro- 
duce relativistic jets. In order to account for the ob- 
served emission, we only need an initial Lorentz factor of 
Tq > 2. Such jets would generally not produce a GRB, 
which typically requires Tq > 100. In this case, if we 
use a conserv ative estimate , comb ining the 33 SNe from 
the sample of iBerger et al.l l|2003fl . and the additional 7 
(inclu ding 2001em) from the sample of iSoderberg et all 
( 2004), then SN 2001em would be one out of 40 nearby 
SNe Ib/c that produce d relativistic jets . Thi s implies 
/rj > 2.5%. Following ISoderberg et all |200l and us- 
ing only nearby SNe Ib/c for which there are late time 
(> 100 days) observations, we obtain /rj ~ w 6.7%. 
Since the observations are sparse (and in most cases con- 
sist of a single upper limit) the actual value of /rj might 
even be larger. 

It is interesting to compare /rj to the fraction /grb 
of SNe Ib/c that produce GRBs. There are various es- 
timate s for /orr . Assu ming a uniform jet with sharp 
edges, iFrail et all l)2001l) found a beaming correction of 
(/ b _1 ) ~ 500 between the observed and the true GRB 
rates (where h ~ fln/2) that results in /grb ~ 0.4%. 
iPerna eta l. (2003) estimated /grb for the universal 
structured jet (US J) model and found /grb ~ 8 x 10~ 6 . 
iGuetta et all (|2003|i found that the USJ model is not 
consistent with the observed log N — log S distribution, 
and did a more thorough analysis for the uniform jet 
model, that resulted in (f^ 1 ) ~ 75 ± 25 and /grb ~ 
(5.5 ± 1.8) x 10~ 4 . Therefore, if indeed / RJ > a few 
percent, then /rj//grb ~ 10 2 , implying that SNe Ib/c 
produce ~ 100 times more mildly relativistic jets (with 
To > 2) than highly relativistic ones (with To > 100) , 
as suggested by seve ral a uthors ||M^F^lveii^t_a l.ll2001t 
lR.amirez-Ruiz et al.ll200l IGranot fc Loebll2003fl . 

If this is the case, one might expect a smooth and con- 
tinuous distribution P(Tq) of initial Lorentz factors for 
the jets produced by SN Ib/c, where Tq > 100 would pro- 
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duce a GRB, Tq > 10—20 could result in X-ray orphan af- 
terglows and possibly also X-ray flashes, To > 5 — 10 may 
give rise to optical orphan afterglows, and To > 2 could 
be responsible for radio orphan afterglows. If, for ex- 
ample, we parameterize this distribution as a power law, 
P(r ) = KT Q v for r min < T < r max , where T min w 1, 
T max > 100, and 4 K = / RJ (r? - 1)1^, then we need 
j] ~ 2 in order to get /rj//grb ~ 10 2 - However, this is 
still highly speculative at this stage. 

We now comp are /rj//grb ~ 10 2 to current obser- 
vational limits. iLevinsion et alJ l)2002l) have used the 
FIRST and NVSS surveys to place limits on orphan radio 
GRB afterglows. They estimated the number of candi- 
dates for such events over the whole sky above 6 mJy to 
be 227, and obtained a lower limit on f^ 1 ~ E~ i so /E, 
of ) > 13. However, this was derived assuming 
a fixed value for the isotropic equivalent energy out- 
put in 7-rays, -E 7l i SO , while allowing the true energy 
E to vary. If instead we fix the true energy t o be 
E k, 10 5 erg, as sugg ested by iFrail et all 0)01) and 
iBloom et alJ ()2003fl , the same analysis would result in 
an upper limit of (J" 1 ) < 6300. Following iGuetta et all 
lt200l . /rj//grb ~ 10 2 and {f^ 1 } = 75 ± 25, which 
is consistent with the revised limit of (/ b _1 ) < 63 that is 
obtained by scaling up the expected number of such tran- 
sients by a factor of /rj//grb- This g ives roughly the 
right n umber of radio transients found bv lLevinsion et all 
( 2002), if most of them are caused by T > 2 jets pro- 
d uced in SNe Ib/ c. 

iNakar &: PirarJ l)2003|) estimated the ratio of on-axis 
orphan X-ray afterglows (T > 10 - 20) and GRBs (r > 
100) to be less than 8, using the ROSAT all sky survey. 
This is marginally consistent with r\ ~ 2 and suggests rj < 
2. Finally, we note that even if the radio emission from 

4 Here J P(To)dTo = /rj is normalized to the total fraction of 
SNe Ib/c that produce relativistic jets. 



SN 2001em arises from the deceleration of a relativistic 
jet, then there is still a large statistical uncertainty on 
the value of /rj, since it is estimated on the basis of one 
event. For example, /rj//grb might still be ~ 10, which 
would imply rj ~ 1.5. 

A relatively large value of /rj might support the 
idea that at least some core collapse SNe, and in par- 
ticular SNe Ib/c, m ay be triggered by bipolar jets 
(|Khokhlov et alJll99flt> . Even if only - 1% of the core 
energy is channelled into such jets, they would still have 
enough kinetic energy to provide most of the power in 
the explosion, and substantially alter the structure of 
the expanding SN shell. While most rotating magnetized 
proto-neutron stars with low power are expected to pro- 
duce broad slowly collimating jets, a few high pow er ones 
shoul d produce n arrow rapidly collimating jets ljUsov I 
H991 lThomp^oiTIIT99l . Although carrying more power, 
these highly collimated jets will be much less efficient 
than the broad jets i n imparting energy an d momentum 
to the outer layers (Khokh Tov et al.l 11999ft. They may 
then act similarly to th e failed SNe l)MacFadven et alJ 
120011 llzzard et alll2004|) . continuing to accrete much of 
the surrounding stellar layers and collapse to a black 
hole, potentially resulting in even faster and narrower 
jets. Observational estimates of the ratio /rj//grb will 
be valuable for constraining the different stellar evolution 
routes involved in producing relativistic, bipolar jets in 
core collapse SNe. 
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